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Population and phase dynamics of F = 1 spinor condensates 
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We show that the classical dynamics underlying the mean-field description of homogeneous mix- 
tures of spinor F = 1 Bose-Einstein condensates in an external magnetic field is integrable as a 
consequence of number conservation and axial symmetry in spin space. The population dynamics 
depends only on the quadratic term of the Zeeman energy and on the strength of the spin-dependent 
term of the atom-atom interaction. We determine the equilibrium populations as function of the 
ratio of these two quantities and the miscibility of the hyperfine components in the ground state 
spinors are thoroughly discussed. Outside the equilibrium, the populations are always a periodic 
function of time where the periodic motion can be a libration or a rotation. Our studies also indicates 
the absence of metastability. 

PACS numbers: 05.30.Jp, 03.75.Hh 
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I. INTRODUCTION 



II. THE EQUATIONS OF MOTION 
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The recently realized trapping of sodium atoms by 
purely optical means opens up the possibility of study- 
ing "spinor" condensates in which the spin degrees of 
freedom are not frozen 0. 

Many authors have investigated, in the framework of 
Bogoliubov theory, the ground state configurations and 
the low-lying collective excitations of homogeneous mix- 
tures of = 1 spinor condensates in the absence 3] and 
in the presence of an external magnetic field 2, 4J . These 
studies predicted a variety of new phenomena such as the 
existence of spin domains in the ground state 0, Q and 
the propagation of spin waves 0, |j| . 

The mean-field dynamics inherent in these works is 
known to be equivalent to a classical dynamics whose 
degrees of freedom are the phase and population of the 
three hyperfine levels 3, 4]. 

In this paper we show that this classical dynamics is 
integrable as a consequence of number conservation and 
axial symmetry in spin space. By a proper canonical 
transformation it reduces to a dynamics involving only 
one degree of freedom. The determination of the equi- 
librium points reveals a rich structure in phase-space. 
Contour plots of the constant energy surfaces show that 
the population is a periodic function of time, where the 
periodic motion can be a libration or a rotation. 

Our studies differs from references II 3 by taking 
explicitly into account the constraint of the axial sym- 
metry in spin space. This allow us to make a detailed 
discussion of the properties of the equilibrium configura- 
tions and of the population and phase dynamics of the 
spinor condensate, which complements previous works 



The hamiltonian of our system of = 1 homogeneous 
mixture of bosonic atoms in the presence of an external 
magnetic field is equal to 0, S 01 
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In |(2Jl, - creates an atom in the hyperfine level a, 

a = 1, 0, —1, with momentum hk, p and q are the inten- 
sities of the linear and quadratic terms of the Zeeman 
energy and cq and C2 are, respectively, the strengths 
of the spin-independent and spin-dependent terms of the 
atom-atom interaction 0, H, Q . 

The hamiltonian (j^J is number conserving and axially 
symmetric in spin space. 
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In the mean-field theory we suppose that the conden- 
sate is a coherent combination of atoms in the p — 
state, 



0) 
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where |0} is the vacuum. The complex numbers are 
the condensate wave-functions for the atoms in the hy- 
perfine level a. To find the time evolution of Zq,, we use 
the time-dependent variational principle 0, |^ 
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Two of these variables are, respectively, the mean num- 
ber of atoms 



SS I ih{z\z) - {z\H\z)dt = 
which reduces to 
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Ni = (z|7V|i) 
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and one-half the mean value of the component of the total 
hyperfine spin in the direction of the magnetic field 



= S J ihz*^Za - Ho{z, z*)^dt = Q (4) 
where the hamiltonian Hq{z, z*) is given by 
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In terms of these new canonical variables the action 
becomes equal to 



Ho{z,z*) = {z\H\z) = -p {\zi\'^ - \z^i\ 



+ 2|zop(i^ip + k-in 

+ 2zi*zliZo2 + 2ziz_i(z*)2) 
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Imposing that the action is stationary with respect to 
variations of Zq, we get Hamilton equations of motion in 
complex coordinates 



ihZa 
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To take advantage of number conservation and axial 
symmetry in spin space we perform two canonical trans- 
formations. The first one introduces, as canonical vari- 
ables in phase-space, the population and phase of each 
hyperfine component by the transformation 



such that the action reduces to 



^ = (^ft^ ^a^a - Ho{N, dt (8) 

The second canonical transformation is given by 
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E HnJ^ - Ho{N,,N2, iVa, §2) dt (12) 



where the Hamiltonian is ciclic in the coordinates 0i and 
^3, 
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As stated before, the dynamics follows from the condi- 
tion that the action (|12|l is stationary, which leads to 
Hamilton equations of motion in these new canonical 
variables, 



dHo 
dNi 
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The property that the hamiltonian is ciclic in 9i and 0^ 
has two important consequences. One is that the mean- 
value of the total number of particles and of the compo- 
nent of the total spin of the condensate in the direction of 
the magnetic field are constants of the motion, Ni = N, 
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2tiN3 = {z\Sz\z). The other is that the dynamics involves 
only one degree of freedom 



(15) 



To simplify the equations of motion (|15|l we define new 
variables which are the old variables divided by the num- 
ber of particles, fia = to write the equations of mo- 
tion as 
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where p = y is the density of the condensate. 

The equations Hlt)|) are the analog of the mean-field 
classical equations of motion that describes atoms in two 
states coupled by a Josephson-type term 0- 

From equations (|16|l . the following general properties 
of the equations of motion emerges: 

1) The population dynamics is independent of the 
strength of the linear term of the Zeeman energy and 
of the spin-independent component of the atom-atom in- 
teraction. Indeed, only the phase 63 depends on p and 
only the phase 9i depends on cq. 

2) By a proper choice of time scale, r = \itp\^ 
population dynamics, in the limit g = 0, is independent 
of the magnitude of C2, depending only on its sign. In 
general, that is when g ^ 0, the population dynamics 
depends on the ratio ^ and on the sign of C2- 



III. PROPERTIES OF THE EQUATIONS OF 
MOTION 



A. Equilibrium configurations 



In equations 1)16(1 and (|17all7b|l the constant of the 
motion n^, is defined in the interval — i < < j and 
the dynamic variable n2 in the interval — ^ < 71-2 < | — 
2 1 77,3 1 . In our discussion of the solutions of the equilibrium 
equations we consider separately the cases ^ Q and 

^3 = 0. 

a) n3 ^ 

We have solutions which depends on the phase 62 and 
solutions which are independent of 62 

al) Solution which depends on O2 with cos 202 = 1- 

In this case the equilibrium value of n2 is given by 
equation (|17b|l with cos 2^2 = 1 • This equation have one 
solution in the interval — 00 < ^ < 1 + \J\ — {2n^Y. 
When ^ — > — 00 the equilibrium value of n2 approaches 
the upper boundary, n2 — |— 2|n3|, for which the fraction 
of atoms occupying the hyperfine levels are ni — |n3|-|-n3, 
Tip = 1 — 2|f7,3|, n_i = |n3| — 713. On the other hand when 
^ = 1 + ^Jl — {2nzY, ^2 is at the lower boundary, 
n2 = —\^ in which case the fraction of atoms occupying 
the hyperfine levels are ni — ^(1 + 2n3), — 0, ri-i — 
\{l-2n^). 

When we neglect the quadratic term of the Zeeman 
energy, that is q = 0, the equilibrium value of n2 is, 
77-2 = \ {\ — (2n3)^) and the occupation fractions are 

2 ._M|i3)!^,_^= i(i_2n3)^ 



ni 



= i(l + 2^3)2, 7io 



a2) Solution which depends on 62 with cos 26*2 = — 1 

The equilibrium value of 77.2 is now given by equation 
(|17b|l with cos 2^2 = ^ 1 • This equation has one solution 
in the interval -2- > 1 - Jl - {2n^Y . When = 

C2P V \ 0/ g^p 

1 — \Jl — (277.3)^ 77.2 is at the lower boundary 77,2 = — ^. 
On the other hand when ^ cxj it approaches the 

upper boundary 712 = | — 2|77,3|. 



From the equations of motion, equations (|16|l . we see 
that the equations which determine the equilibrium con- 
figurations are 



3-7.2 
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a3) Solution which does not depend on the phase O2 



In this case 77,2 is at the lower boundary 772 = —3, 
independently of the value of 

These properties are illustrated in FIG^lwhere we plot 
the equilibrium values of 77,2 as a function of ^ for 773 — 
1 

I' 



b) 773 = 

In this case the equilibrium equations, H17all7b|l . re- 
duces to 
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TABLE I: Equilibrium configurations for different parameter domains, 7^ 0, antiferromagnetic case, C2p > 0. 
< 1 — \/l — (2ri3)^ al) maximum; a3) minimum 

1 — \/l — (2n3)^ < j^^^ < 1 + \/l — {2n:i)^ al) maximum; a2) minimum; a3) undefined" 

> 1 + \/l — (277,3)-^ a2) minimum; a3) maximum 

"The classification undefined means that when we leave the line 
defining the corresponding boundary the energy increases or de- 
creases depending of the value of the phase 52- 

TABLE II: Equilibrium configurations for different parameter domains. 723 7^ 0, ferromagnetic case, C2p < 0. 

al) minimum; a3) maximum 

al) minimum; a2) maximum; a3) undefined 

a2) maximum; a3) minimum 



- (1 + vT^T2^) < ^ < - (1 - \/l - (2^3)^) 
^<-(l + Vl-(2n3)^) 



n2 




FIG. 1: Equilibrium values of n2 as function of -2- for 

^ C2P 

— J. Solid curve is al), dashed is a2), and the lower 
boundary n2 = — | is a3). The horizontal straight lines are 
the boundaries of n2. 
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C2P \ 3 
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sin 26*2 = 
OS 2^2) = 



Again, we have two classes of solutions, dependent and 
independent of the phase 62 ■ 

bl) Solution which depends on the phase 62 with 
cos 2^2 = 1 

The equilibrium value of n2 is n2 = ^ — 2^) 
interval —2 < ^ < 2. In this case the occupation frac- 
tion of the hyperfine levels are ni = n_i = i (l + 



no = 



2 2C2P J 



b2) Solution which depends on the phase 62 with 
cos 202 = -1 



In this case there is a solution only at g = and it is 
the hue cos202 = —1, — < ^2 < |. 

When 713 = 0, there are two solutions which does not 
depend on the phases. 



b3) One is the lower boundary 712 = — | in which case 
the occupation fractions are ni = 77_i = i no = 0- 



b4) The other is the upper boundary 77,2 = | in which 
case the occupation fractions are Tii = n_i — 0, no — 1. 
The solutions b3) and b4) exist for any value of 

The properties of the 77,3 = equilibrium configurations 
are illustrated in FIG|21 , where we plot the equilibrium 
values of fia as a function of -2-. 

We would like to point out that it can be shown that 
the phase dependent solutions of the equilibrium equa- 
tions are the roots of the third order equation. 



2^ /0 

C2PJ 



4 ^/o^--^ 2 + ^ /0^ -(2773)^2 1 

C2P C2P \ C2P ; 



{^n,f 1 + ^ 

\ V C2P 



(2773) 



/o7^2|773|(18) 



with 772 = f ~/oj where /o is the fraction of atoms outside 
the 777/ = hyperfine level, Jo = n\^ n^\. 

We can find analytical expressions for the roots of 
equation H18|) but they are not particularly illuminating 
and we will not write them here. 
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TABLE III: Equilibrium configurations for different parameter domains, na — 0, antiferromagiietic limit, C2p > 0. 



< -2 



-2 < < 



< < 2 



r^>2 



b3) minimum; b4) maximum 
bl) maximum; b3) minimum; b4) undefined 
bl) maximum; b2),b3),b4) degenerate minimum 
bl) maximum; bS) undefined; b4) minimum 

b3) maximum; b4) minimum 



TABLE IV: Equilibrium configurations for different parameter domains, ns = 0, ferromagnetic limit, C2p < 0. 



< < 2 



g = 



-2 < < 



< -2 



b3) maximum; b4) minimum 
bl) minimum; b3) maximum; b4) undefined 
bl) minimum; b2),b3),b4) degenerate maximum 
bl) minimum; b3) undefined;b4) maximum 
b3) minimum; b4) maximum 



n2 
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\ 

\ 

0>2 
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FIG. 2: Equilibrium values of n2 as function of ^p- for na — 0. 
The dashed straight line is bl). The vertical straight line 
Q = is b2). The horizontal straight lines are the boundaries 
of n2. 



A summary of our discussion is displayed in the tables 
(|IIIIIIIIIIVI| ■ where we consider separately the antiferro- 
magnetic and the ferromagnetic limits. From the tables 
we can easily determine what are the ground state config- 
urations for the different parameter domains and in the 
two limits considered above. 

Ground state configurations have been investigated in 
reference 0, in the absence of the magnetic field {p — 
q — 0). They observe that the (degenerate) ground state 
in the antiferromagnetic limit are the polar states and in 
the ferromagnetic limit, the ferromagnetic states. 

The polar states of Q have 77.3 = and are the states 



b2),b3), and b4) at g = (see Table III). When = 0, 
the ground state in the ferromagnetic limit is the state 
bl) which is equal to the ferromagnetic state of Q with 
n3 = 0. 

When ^ 0, the states al) at g = are equal to 
the ferromagnetic states of (3| and are the ground states 
in the ferromagnetic limit (see Table II). On the other 
hand, the ground state in the antiferromagnetic limit is 
the state a3), which is not a polar state. 

Reference investigate the ground state configura- 
tions in the presence of the magnetic field and neglecting 
the quadratic term of the Zeeman energy. As shown in 
this paper, in this case the ground state configurations 
coincide with the ground state configurations in the ab- 
sence of the magnetic field, a fact overlooked in Be- 
sides a state with maximum value of 713, = ^, 
identify only the spinor a3), which is the ground state 
just in the antiferromagnetic limit (see Table I). 

The general case when we consider both the linear and 
quadratic terms of the Zeeman energy have been stud- 
ied in reference . Our approach differs from 2] , in the 
sense that we take explicitly into account the constraints 
of axial symmetry which simplifies considerably the dis- 



B. Miscibility 

One question that we can adress is the miscibility of 
the hyperfine components in the ground state spinors . 

Consider first the case ^ and the antiferromag- 
netic limit. For q < I — y/l — (277.3)2, the ground state 
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FIG. 3: Top: energy surface as function of 62 and n2- Bot- 
tom: contour plot of the energy surface in the 62 x n2 plane. 
In these plots na = i and ~ 0.4, in the antiferromag- 

netic limit. The energy is in units of \c2p\. Darker colors 
means lower energy. 



spinor is a3), where only the m/ = ±1 components are 
miscible , the population of the my = component be- 
ing null. For > 1 — \Jl — {2nzY , the ground state 
spinor is a2), in which case the three components are 
generally miscible. Actually the miscibility of the hy- 
perfine components depends on the ratio f^^- When 
approaches its lowest value, n2 is near the lower 

boundary, 712 = —\, therefore the rrif — Q component 
is practically immiscible with the ruf = ±1 components. 
When increases, a2) approaches the upper bound- 
ary n2 = I — 2|n3|, in which case the m/ = component 
mixes predominantly with the nif = 1) component 
depending on the sign of 713, > 0(< 0), the population 
of the third component being negligible. 

The next case is fia 7^ and the ferromagnetic 

limit. For < ^ (l + v^l - (Zns)^^ , the ground 

state spinor is a3), where only the rrif = ±1 compo- 
nents are miscible, the population of the mf — com- 
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FIG. 4: Top: energy surface as function of 62 and n2. Bot- 
tom: contour plot of the energy surface in the 62 x 712 plane. 
In this plots ns = and = 0, in the ferromagnetic limit. 
The energy is in units of \c2p\. 



ponent being null. On the other hand, when > 

— ^1 + \/l — (277,3)^^ the gromid state spinor is al). As 

in the antifcrromagnetic limit, in this case the three hy- 
perfine components are generally miscible, the degree of 
miscibility depending on the ratio j^^- When this ratio 
approches its lowest value, n2 is near the lower bound- 
ary, n2 = and the nif = component is practically 
immiscible with the = ±1 component. On the other 
hand, when -r-^^ increases approachig -l-oo, n2 is near the 

upper boundary, 712 = | — 277,3, in which case the mj = 
component mixes predominantly with the mf ^ +1(— 1) 
component depending on the sign of 77,3, 773 > 0(< 0), the 
population of the third component being negligible. 

The only case left is when 77,3 — 0. Consider first the 
antifcrromagnetic limit. For < the ground state 

O \C2P\ ° 

spinor is b3), for which 712 is at the lower boundary, 
77,2 = — -I . In this spinor we have equal population of 
the ruf = ±1 components, the population of the mf = 
component being null. On the other hand, for j-r-r > 0, 

I C2p 
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712 ^2 




(c) (d) 

FIG. 5: (a)n2 ,(b)^2 as function of time, in the case of rota- 
tion; (c)n2, (d)^2 as funtion of time, in the case of libration. 
Notice that 712 is always periodic. The time is in units of 
and the angles, in radians. The parameters are equal to the 
ones in FIG|3| The initial conditions are: n2 = 0.159 and 
62 = J for rotation, ?i2 ~ —0.1 and ^2 = f for hbration. 

the ground state spinor is b4) for which n2 is at the 
upper boundary, n2 = |. In this spinor all atoms are in 
the nif = state. 

The ferromagnetic limit is richer than the previous one. 
For < —2 the ground state spinor is b3). How- 

ever, for —2 < < 2, the ground state spinor is bl), 
where the three components are generally miscible. Ac- 
tually, it changes from a spinor where practically only 
the my = ±1 are miscible, the population of the ruf — 
component being negligible, near the lowest value of -j^^ 
to one where practically only the component my = is 
populated near the highest value of For > 2 

the ground state spinor is b4). 

IV. DYNAMICS 

The qualitative features of the population dynamics 
can be easily visualized if we make portraits of the 



contour curves =constant in the phase space plane 
92 X n2- Examples are shown in FIG|21and FIG0] We 
see that n2{t) is always a periodic function of time. The 
motion can be a libration when 62 is a limited function 
of time and a rotation in which case 62 always increases 
(decreases) as a function of time. Examples of these be- 
haviors are displayed in FIGEI 

V. SUMMARY AND CONCLUSIONS 



In summary we have studied in this paper the classical 
dynamics that underlies the mean-field description of an 
homogeneous mixture of spinor F = 1 condensate in an 
external magnetic field. 

As a consequence of number conservation and axial 
symmetry in spin space this dynamics is integrable. The 
equations of motion show that the population dynamics 
depends only on the quadratic term of the Zeeman energy 
and on the strength of the spin-dependent component of 
the atom-atom interaction. 

For a fixed mean-value of the component of the con- 
densate spin in the direction of the magnetic field we 
determine the equilibrium configurations as a function 
of the ratio -2-. We also make a detailed discussion of 

C2P 

the miscibility of the three hyperfine components in the 
ground state spinor as a function of Our studies re- 
veal the absence of metastability in the sense that there is 
no two local minima configurations in the parameter do- 
main. We have shown that outside the equilibrium, the 
populations are always a periodic function of time, where 
the periodic motion can be a libration or a rotation. In 
the first case the phase 62 is always limited whereas in 
the second case it always increase (decrease), as function 
of time. 

Finally we would like to remark that the restriction to 
a homogeneous mixture is not only of academic interest 
0, 0]. Besides being a guide to what happens in the 
case of trapped condensates, classical Hamiltonians of 
the type considered in this paper emerges in a mean-field 
description of the quantum single-mode approximation 
for spinor condensates in a trap 8| . 
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